Slave fermion approach to the metalHc phase in large Ud high-Tc cuprates 
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Atomic repulsion Ud on the Cu site in high Tc cuprates is large but, surprisingly, some important 
properties are consistent with moderate couplings. The time dependent perturbation theory with 
slave particles is therefore formulated in the Ud ^ oo limit for the metallic phase in the physically 
relevant regime of the three-band Emery model. The basic theory possesses the local gauge invari- 
ance asymptotically but its convergence is fast when the average occupation of the Cu-site is small. 
The leading orders exhibit the band narrowing and the dynamic Cu/02 charge transfer disorder. 
The effective local repulsion between particles on oxygen sites is shown to be moderate in the phys- 
ical regime under consideration. It enhances the coherent incommensurate SDW correlations. The 
latter compete with the Cu/02 charge transfer disorder, in agreement with basic observations in 
high Tc cuprates. 
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It has long been maintainedii^ and recently 
reemphasized^ that one of the most important 
unanswered questions in the theory of the high-Tc 
superconductivity in cuprates concerns the nature of 
metallic correlations which reduce the large Coulomb 
interaction Ud on the Cu sites. In the present work we 
discuss this issue within the widely used large Ud Emery 
model^ which adds two oxygen sites Ox,y to each Cu 
in the Cu02 unit cell. The large Ud limit is usually 
taken perturbatively in the Cu-0 hybridization, starting 
from the unperturbed state with n^^'' = 1. This leads 
to the superexchange^ J as the basis of the popular 
tJ models^. However, various local measurements for 
doping < 0.2, NQR in particular, indicate^ that the 
average occupation of the Cu site is n^j « 1/2. We argue 
here that the nj^"'' = unperturbed ground state is a 
good starting point for the large Ud perturbation theory 
associated with the metaUic phase at appreciable |a;|, 
while the x ~ Mott-AF phase and the short range AF 
order can be approached from n^"-* = 1 (tJ model) side. 
Thus we generaHze to finite Cu-0 hybridization the early 
proposal^ that in the high Tc cuprates the d-level with 
large Ud falls close to the nearly free (oxygen) band and 
that it can be related to the dio dg disorder effects. 
In particular we show that, in the presence of the Cu-0 
hybridization, the dio dg disorder is dynamic, while 
the effective interaction between the carriers propagating 
on oxygens is repulsive and reasonably small. 

As usuati, the Cu,0 site energies are denoted here by 
p — Ed > in the hole picture. The 
-Oy hoppings are tpd and tpp, the latter 
being retained although smaller than tpd. Finally 1 -I- a; 
denotes the total number of holes per Cu02 unit cell. 
The equal sharing of the hole charge between Cu and 
two 's in the metallic phase is characteristic of the free 
fermion tpd ^ Ap^ model^"^. The values « 1/2 were 
also shownifiiii to be consistent with the mean-field slave 
boson (MFSB) fits of the ARPES results, which suggest 



This motivates us to investi- 
Emery model 



Ed, Ep and Apd 
Cu-0 and 



the regime 2i^^ < -Apdtpp 
gate here the Ud — oo limit of the finite 
for Apd sufficiently small and tpd > \tpp\. 

To this end, we use the Ud = oo slave particle 
theoryiSiii^'i^iii. The three-state space {d^ state omit- 
ted) associated with the i?-th Cu site is described using 
the f \ and 6t slave particles and the physical fermion 

ct, is represented as ct; &t; fs- The anticommu- 

tation rules of the physical fermions ct on and among 

/t,cr 

Cu sites are satisfied, provided that the /t and 6^ ^ slave 

particles are taken, respectively, as fermions and bosonsi^ 
(slave fermion theory; SFT) or vice versaii^ (slave boson 
theory; SBT). In both representations the number op- 
erators satisfy Qr ^ nf^-^ + X^o- ^f>H,cr = 1. The corre- 
sponding number operators obey ndji,a = nb^,a-, which 
is usually calledi^ the Luttinger sum rule (LSR). The 
anticommutations between the Cu and 0-sites are sat- 
isfied in the SBT because 6^ and correspond to 

indistinguishable fermions. However, in the SFT they 
are replaced by commutations because this theory deals 
with three kinds of distinguishable particles. Therefore 
the SFT has to be antisymmetrized a posteriori. 

The effective repulsion in the Ud ^ oo limit is a ki- 
netic, time-lag effect^ i^^i^^ in the sense that one particle 
has to wait for the other one to leave the Cu site in order 
to cross it. It is therefore appropriate to use the time- 
dependent theory. The convenient operative choice is the 
time-dependent many-body perturbation theory which 
treats tpd to the infinite order of perturbation. In this 
theory it is essentiaU^ to start from the nondegenerate 
unperturbed ground state of the tpd ~ Hamiltonian 
Hq. The slave particle representation via slave bosons or 
slave fermions has to be chosen accordingly and we turn 
first to this question. 

The unperturbed slave particle ground state |Go) has 
itself to be metallic and at tpd = it is the direct product 
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of the states of three kinds of particles, \gI)(E)\Gq)(E)\Gq) . 
The only way to ensure translational and local gauge in- 
variance in the unperturbed state is by taking = 0. 
|Gg) is then the usual Hartree-Fock (HF) state of 2?!^°^ = 
1 + a; fermions p-''^ in the cosine band associated with 
Ep and tpp, which parametrize Hop. This band is filled 
up to the chemical potential f^i+x and is usually folded 
artificially in two i ~ 1,1 oxygen bands in the Cu02 
Brillouin zone. Such a |Gq) is, of course, nondegenerate 
in both slave particle representations. The n^p — 1 d^" 
state on Cu is denoted by /j^jO) , where |0) is the aux- 
iliary vacuum on Cu. IGg) is a nondegenerate state of 
the tpd ~ slave particle Hamiltonian Hqx only if fj^- 
particles are taken as spinless fermions rather than as 
spinless bosons. Only then can |Gq) be equally simply 
expressed in terms of Fourier transforms fl of the oper- 
ators /Jj since it is only then that we have 

iG^)=n4io>=n4io>' (1) 

R k 

up to an unimportant phase factor. Here the product 
over k covers the whole Cu02 Brillouin zone. In other 
words, the Mott state of spinless fermions is equivalent to 
the filled (dispersionless) band of those particles. When 
the ft are chosen as spinless fermions, |Gq) is the state 

with no &-bosons, i.e. n^°^ = or, equivalently, the (P 
states fet^ 1 6) are absent in I Go). The overall unper- 

turbed n^°^ = n^^^ = translationally and locally gauge 
invariant ground state is thus nondegenerate. The price 
of this important achievement is, however, that the SFT 
must be antisymmetrized a posteriori. 

Equation (U) represents the key step of the present ap- 
proach. Once the unperturbed ground state |Go) is ex- 
pressed in the momentum space so too should the transla- 
tionally and locally gauge invariant slave particle Hamil- 
tonian Hx = Ha - XN + Hj, Ha = H^p + Hqx, 

i,k,(7 k,CT k 

1.(7, k.q 

t^ik) = t,,y2(^|sin^|±|sin^|) . (2) 

where N is the number of the Cu02 unit cells. The 
actionii of Hj on the state of Eq. Ilj is then simply 
defined in the full momentum representation. Provided 
only that it is convergent, the ensuing time-dependent 
slave fermion perturbation theory (SFT), in terms of 



Hj, generates then asymptotically the exact SFT ground 
state, which is translationally, time inversion and^^ lo- 
cally gauge invariant. We recall however that such a 
procedure replaces the Cu-0 anticommutations by com- 
mutations. 

This omission is, however, irrelevant in the lowest or- 
der of the Dyson perturbation theorjJi. The reason is 
that the Cu site is initially empty, while one Cu hole is 
required for the anticommutator (and two for Ud inter- 
action). The r = 1 pdp or dpd propagators (particles 
created and annihilated on O or Cu sites respectively) 
correspond to the lowest order irreducible Dyson p-self- 
energy proportional to the convolution B^^'' * pf'' of the 
bare slave particle propagators. If the chemical poten- 
tial of the p-fermions pS^^ is determined subsequently 
from the relation ri^p + 2rL"p^ = 1 -I- a;, keeping in mind 
that ultimately Ud = nt in the SFT, the pdp and dpd 
propagators become equivalent to those of the tpd hy- 
bridized HF theory. The HF result is characterized by 
three i = L,I,U branches of poles to), of the free three- 

k 

bandit model, with large Fermi arcs associated with the 
lowest L-band. The corresponding spectral weights zl' 

are given in terms of the departure of , i = L,I, U, 

from the unhybridized Ed and ei'\ i = 1,1. For exam- 
ple, the spectral weight of the L-band in the r = 1 pdp 
propagator is 

, , (4^)_£,)2(4^)_,(0)(4^)_/)) 

(L) _^ k °-' ^ k pk'^ k pk' 

^Pdy^k ~^k '^'^k ~'^fe ) 

The normalization of zl is chosen so that the hop- 
ping from the 0-sites to the Cu-site is associated with 
tpd, i.e. the fc-dependences of tpj(fc) are absorbed in z'j^K 

It is noteworthy that, due to anticrossingi^, the zl*-* are 

regular functions of k. The lowest order SFT generates 
the HF propagators irrespectively of the average HF oc- 
cupation n^^^ of the Cu-site. On the other hand, it will 
be seen below that n'"^^ w n^^'' only when n^J'' is small 
and similarly then Q^^) = nf^ + n^^^ w 1. The SFT thus 
converges quickly only when n^l^ is small. It is therefore 

important to keep in mind thalii n|;^^ w 1/2 at |a;| small 
for tpd ^ Apd and that finite tpp decreasesi^ it further. 

The useful overall rule of thumb is that n''p < 1/2 as 
long as the HF chemical potential /i'-^-' falls below the 
vH singularity at lOvh in the lowest L-band i.e., as long 
as a; < XvH wherein x^h ~ —tpp is the positive doping 
required to reach the vH singularity. 

The central quantity for further expansion is the 
effective interaction between the p-particles shown in 
Figs. [lK,b. It expresses the fact that two holes can- 
not hop simultaneously to the same Cu-site. When 
the ipd-anticrossing of the d-level and the oxygen band 
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Figure 1: Effective scattering vertex of two pdp propaga- 
tors. Triangles denote multiplication by tpd. Internal blue 
and green lines are the dispersionless free propagators of the 
b- and /-particles in the state of Eq. |[T|, respectively, with 
appropriate arrows of time. 



strongly affects the occupied states in the p-bands, it 
is appropriate^^, as assumed hereafter, to consider the 
processes in Fig. [T] as the interactions between the r = 
1 hybridized pdp-propagators normalized according to 
Eq. |[3|. The triangular vertices in Figs. [T] are then to 
be associated with tpd- The internal arrows of time ex- 
plicitly account for the fact that the bosons can only be 
created and the spinless fermions only annihilated in the 
unperturbed ground state in Eq. ([TJ . This defines unam- 
biguously the temporal structure of the internal squares 
AJ5 and Ao in Figs[Ih,b. E.g. Ajj of Fig.[ljD is given by 



A 



2ed -uji - UJ2 



ill) ' 



(4) 



Here ujj denote the external frequencies which run coun- 
terclockwise around Aq in Fig. [TJd with ui in the upper 
left corner. The result ^ is independent of A in Eq. |[2]). 
This illustrates one general satisfactory featurei^ of the 
perturbative SFT, namely that it is independent of A or- 
der by order even though the local gauge invariance is 
only asymptotical. 

It is important to note, however, that spin is conserved 
in the triangular vertices of Figs. [lh,b i.e., that in the 
SFT two pdp-particles cannot hop simultaneously to the 
same Cu site through the assistance of slave particles, 
irrespective of their spin. Contrary to that, for the origi- 
nal Emery Hamiltonian with Ud, it is the Pauli principle 
which forbids two particles in the triplet configuration to 
occur simultaneously on the Cu-site irrespective of Ud- 
The action of the local two-particle interaction on the 
fermions of the same spin vanishes thus identically after 
antisymmetrizationi^. This discrepancy with the SFT 
stems clearly from the omission of the Cu-0 anticom- 
mutations in the SFT and the triplet scattering has to 
be remove d^^i^^ from the SFT in order to make it cor- 
respond to the original problem. The SFT modified in 
this way, i.e., antisymmetrized a posteriori, is named the 
MSFT here. 

There are, however, processes, additional to the lowest 
r = 1 HF hybridization, for which this Cu-0 anticom- 
mutation is irrelevant. This is best seen upon closing 
one tpd pdp-hybridized line in each of Figs. [lK,b. The 



effective interaction Aq''' is then absorbed in the Dyson 
bubble renormaHzation of the slave particle propagators 
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Figure 2: Contributions to irreducible p-self energy of the 
r = 4 pdp propagator: (a) in the particle-particle channel; (b) 
in the particle-hole channel. 



the Cu-0 anticommutations. In addition, those propa- 
gators are local but each exhibits temporal disorder—, 
b ^ f symmetrically. As easily seenji^ the associated 



' is small for n 



(1) 

d 



small and it follows that Q'^' 



1. 



The convolution B\' F\' is also local and defines the 
next order (A-independent) irreducible p-self energy. In 
the first place, for n'^ small this gives the band nar- 

= 2 'pd'p and 



rowing 



t'Li^ - "i'V2) in the 



dpd propagations (half of that obtained in the MFSB 
approachesi^^) . In addition the b- and /-disorders con- 
volute (causally) in frequency. This introduces the imagi- 
nary continuum in the r = 2 pdp and dpd propagationsfi^ 
which is causal in time but incoherent in space. In other 
words, this continuum arises from the spatially incoher- 
ent, local, dynamic d^ charge-transfer fiuctua- 
tions. The leading disorder contribution falls in the range 



rrf < < 2/x(^) — Ed, well below the Fermi level. 



S^^'' and F^''', which is insensitive to the omission of 



9, .^^'f 

The clearly distinguishable signature of large Ud in the 
time-dependent perturbation theory is thus the local dy- 
namic c?^" ^ d^ charge-transfer disorder. This can be 
generahzedi^ to r = 3 and further to all processes in 
which the interactions of Figs. [lK,b can be absorbed in 
the bubble renormalizations of the slave particle propaga- 
tors, namely to the non crossing approximation^^ (NCA). 
The broad d^^ ^ d^ continuum is expected then to 
spreads^ all over the spectrum. It is interesting to note 
that a similar broad continuum is obtained in the large 
Ud local2i dynamic mean-field theory (DMFT). 

In contrast, the examples of lowest order, irreducible 
p-self energy terms in which the effective interactions of 
Fig. [T] cannot be dissolved in the bubble renormaliza- 
tions of the slave particles are shown in Fig. [21 With 
other ,f ^ b symmetric terms they define the r = 4 pdp 
and dpd propagators. The SFT and MSFT start to differ 
at this level, i.e. the MSFT requires the omission^^^ of 
the contributions which involve two 6-propagators of the 
same spin within the Aq''' squares. In Figs. [2] this selects 
out the effect of the singlet pdp-pdp particle-particle and 
particle-hole correlations on the single particle propaga- 
tion. Those correlations are best identified on remov- 
ing the upper receding r = 1 pdp-line (together with 
its vertices) from Figs. [21 What remains are the ele- 
mentary pdp particle-particle and particle-hole bubbles. 
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which are then subjected to further renormaHzations by 
tp^Ag'^. The MSFT, amounts thus to the summation 
of a subseries of the locally gauge invariant SFT series, 
and therefore, is not manifestly locally gauge invariant. 
However, if the perturbation procedure for the B^"^^ and 
pir) propagators is carried out b ^ f symmetrically^^ 
the MSFT is at least a conserving^ approximation. 

The SFT tends to the full and the MSFT at least to 
the average^ local gauge invariance in the weak coupling 
manner, provided that ^^^{Aq + A[j) can be treated as 
an effective interaction, small with respect to tpd when 
the external frequencies ujj in Eq. (|4]) fall close to the 

HF Fermi level fi^^K Importantly, small makes the 
effective repulsion t^^Ao small itself. The reason is that 

tpd < £d — fJ-^^^ for fj.^^'' < ujyH (sd — ^vH is linear in 
tpd for smalli Apd and arbitrarj^i^ tpp). According to 
Eq. (0]), the important values of tp^Ag' are then less than 
tpd and, consistently, the internal renormaHzations of the 
slave particle propagators within Aq' can be neglected. 

tpdi^o + ^0) is then the effective interaction between 
the r = 1 pdp propagators. Such a distribution of 
the roles of the propagators and the interactions con- 
trasts with the small Ud theory, where the interaction 
Ud involves the dpd propagators, or with the ansatz 
Ud — > Jpdicoskx + cosky) that is often used2^ in the 
approach from the MFSB side. The r = 1 pdp propa- 
gators with spectral weights of Eq. ([3]), combined with 
tpd{AQ + Aq), resolve thus the long lasting questio n^^i^^ 
of what replaces Jpd ~ Af^^/A^^ in the Ud = 00 metaUic 
limit. 

Further on, the elementary pdp bubble has the same 
nonlocal nesting properties (in the denominator) as the el- 
ementary dpd bubble but differs from it in the form fac- 
tor (numerator The renormaHzations of the elemen- 
tary pdp particle-hole bubble in Fig. [Jb by fj^^^Ao are also 



nonlocal and enhancei^ the coherent SDW correlations 
XSDW- When the elementary pdp particle-hole bubble in 
Fig. [Jb is replaced by such a xsdw, a structure is ob- 
tained which justifies in many respects the semiempirical 
"single loop" approximation used in the calculationii of 
the magnetic pseudogap. In particular the coupHngs of 
the intermittently added pdp carrier to xsdw are thus 
identified here as f^^Ao. 

It can be finally noted that the above choice of the r = 
1 pdp propagators in Fig. [2] is not unique, i.e. that r < 3 
or even the NCA pdp propagators can also be chosen as 
outsets of the renormalization algorithms according to 
the value of x under consideration. The NCA emphasizes 
the c?^" <-*■ c?^ charge-transfer disorder effects while the 
choices r ~ 1,2 favor the coherent SDW correlations. 

In conclusion, tp^(AQ -|-Ag) turns out to be the elemen- 
tary effective local singlet repulsion of two pdp-particles. 
It is the small Apd, n^d^ — counterpart of the nf^ = 1 
superexchange^ Jpd = ^t^^/Ap^ obtained for large Ap^. 
The present perturbation theory is fiexible enoug h^^i^^ to 
deal with the competition between the d^ disorder 

and the coherence associated with the incommensurate^^ 
SDW in simple approximationsi'^'^'^i. Thus the question 
left open here, whether or not the omission of the triplet 
scattering exhausts the a posteriori antisymmetrization 
of the SFT, i.e. whether the asymptotic behavior of the 
MSFT is exactly or approximately locally gauge invari- 
ant, becomes of little practical importance. 
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